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ABSTRACT 

Context. The kHz quasiperiodic oscillations (QPOs) observed in low-mass X-ray neutron star binaries are most likely connected to 
the orbital motion in the accretion disc. The ratio between frequencies of the upper and lower observed QPOs mode cluster close to 
ratios of small natural numbers, most often close to the 3/2 value, but the other rational ratios occur in some sources as well. The 
class of QPOs models considers a resonance between Keplerian and epicyclic frequencies of the geodesic motion. It was suggested 
that a multipeaked ratio distribution may follow from different resonances. The atoll source 4U 1636-53 shows datapoints clustering 
around two distinct values (3/2 and 5/4) of the frequency ratio. The same frequency ratios correspond to the change in the sign of the 
twin peak QPOs amplitude difference, suggesting existence of a resonant energy overflow. 

Aims. We explore the idea that the two clusters of datapoints in 4U 1636-53 result from two different instances of the same orbital 
resonance corresponding to the two resonant points. 

Methods. Assuming the neutron star external spacetime to be described by the Hartle-Thorne metric, we search for a frequency 
relation, matching the two observed datapoints clusters, which may correspond to a resonance. We consider orbital and associated 
epicyclic frequencies with accuracy up to the second order terms in the neutron star angular momentum j and first order terms in its 
quadrupole moment q. 

Results. We have identified a suitable class of frequency relations well fitting the observed data. These models imply for central com- 
pact object in 4U 1636-53 the mass M = 1.6-2.5Mq, dimensionless angular momentum j = 0-0.4, and quadrupole momentum 
q = 0-0.25, with most preferred values M = 1.77M©, j = 0.05, and q = 0.003. 

Conclusions. The relationship implied for a particular case of so called total precession resonance between the Keplerian and the 
total precession frequency introduced in this paper resembles the twin peak QPOs observed in 4U 1636-53 with a^ 2 ~ 3 d.o.f. which 
is about one order lower than x 1 reached by other theoretical relationships which we examine. Moreover, the position of 3/2 and 5/4 
resonant points implied by the total precession relationship well coincides with frequencies given by the change of the rms-amplitude 
difference sign. Notice that if a resonance (in our opinion most likely present) is not considered, the total precession relation has 
a similar kinematic meaning as the periastron precession relation involved in the model of Stella and Vietri, but gives substantially 
better fit and lower neutron star mass. 
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1. Introduction 

Rossi X-ray Timing Explorer (RXTE, Bradt et al., 1993) pro- 
vides observations of the high frequency kilohertz QPOs in the 
X-ray fluxes from neutron-star binary systems (see, van der Klis, 
2006, for a review). 

Several models have been outlined to explain the kHz QPO 
frequencies, and it is mostly preferred that their origin is related 
to the orbital motion near the inner edge of an accretion disc. In 
particular, two ideas based on the strong gravity properties have 
been proposed. While Stella & Vietri (1998, 1999) introduced 
the "Relativistic Precession Model" in which the kHz QPOs rep- 
resent direct manifestation of the modes of a relativistic epicyclic 
motion of blobs in the inner parts of the accretion disc, Kluzniak 
& Abramowicz (2001) proposed models based on resonant in- 
teraction between orbital and/or epicyclic modes related to non- 
linear oscillations of the accretion disc. 

Abramowicz et al. (2003) noticed that ratio between the 
frequencies v L and v v of the lower and upper observed kHz 
QPO mode cluster in neutron-star sources usually close to ra- 
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tios of small natural numbers, most often close to the value 
v u /v L — 3/2. 1 The ratio clustering was later confirmed by 
Belloni et al. (2005). The question whether such a cluster- 
ing represents clear argument supporting a resonance hypoth- 
esis remains to be subject of discussions. Nevertheless it was 
also suggested that, due to multipeaked distribution in the fre- 
quency ratio, more than one resonance could be at work if 
a resonant mechanism is involved in generating the neutron- 
star (Belloni et al., 2005; Torok et al., 2007) and black-hole 
(Stuchlfk & Torok , 2005) kHz QPOs. 

Torok & Barret (2007) realized an interesting root-mean- 
squared-amplitude (rms amplitude) evolution in group of six 
neutron star atoll sources (namely 4U 1636-53, 4U 1608-52, 
4U 1820-30, 4U 1735-44, 4U 1728-34 and 4U 0614+09) - the 
upper and lower QPO amplitudes equal each other when the 
source passes through a rational frequency ratio. Such a be- 
haviour highly suggests the existence of an energy overflow be- 



1 The same 3/2 ratio of QPO frequencies was first noticed in the GRO 
J1655-40 data by Abramowicz & Kluzniak (2001) and later found in all 
other microquasars in which the twin QPOs have been detected (see, 
McClintock & Remillard, 2004). 
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Fig. 1. Left (from Torek et al., 2007): The frequency correlation in the atoll source 4U 1636-53. Curve vk determines the upper QPOs fre- 
quency following from the relativistic precession model (Stella & Vietri, 1999) under the consideration of the gravitational field described by the 
Schwarzschild metric with the central mass M = 1.84M Q . The secondary vertical axes indicates the appropriate dimensionless radius. One should 
note that the assumed 5:4 resonance occurs very close (0.25M) to the innermost stable circular geodesic orbit, i.e., near the expected inner edge 
of the accretion disc. Right (from Torok & Barret, 2007): The rms amplitude difference as a function of the frequency ratio in 4U 1636-53 (black 
points). Gray points correspond to the single detections of the lower (upper) QPOs and the absolute value of their rms amplitude. 



tween the upper and lower frequency mode typical for non-linear 
resonances. The effect seems to be related to the results of the 
frequency ratio distribution studies because the change of the 
rms amplitude difference sign occurs close to the same frequency 
ratio as those corresponding to the datapoints clustering. 

The results of studies mentioned above indicate that for a 
given source the upper and lower QPO frequency can be traced 
through the whole observed range of frequencies but the prob- 
ability to detect both QPOs simultaneously increases when the 
frequency ratio is close to the ratio of small natural numbers 
(namely 3/2, 4/3 and 5/4 in the case of six atoll sources, Torok 
& Barret, 2007). 

Figure 1 illustrates the frequency correlation for the atoll 
source 4U 1636-53 together with its rms amplitude difference 
evolution. 2 For the 3/2 and 5/4 frequency ratio the two dis- 
tinct clusters of datapoints are easy to recognize on the figure 
as well as the corresponding change of the rms amplitude differ- 
ence sign. 

In this paper we focus on 4U 1636-53 and explore the idea 
(Torok et al., 2007) that the two clusters of datapoints mentioned 
above result from two different instances of the same orbital res- 
onance varying in resonant coefficients. 



2. Orbital frequencies of geodesic motion close to 
rotating neutron stars 

In Newtonian physics, a test particle orbits the central mass M 
with the Keplerian angular velocity Q.* = (GM/r 3 ) 1 ^ 2 and its, 
in general eliptic, trajectory is closed, i.e., related azimuthal {or 
Keplerian) v*, radial v*, and vertical v* g orbital frequencies equal 
each other. 

In the general relativistic description, necessary in the case 
of compact objects inluding neutron stars (Misner et al., 1973; 
Shapiro & Teukolsky, 1986; Glendenning, 1997), the trajectory 



2 The data (Barret et al., 2005; Abramowicz et al., 2005) correspond 
to individual continuous segments of observation processed by standard 
shift-add technique (see, Mendez et al., 1998). Notice that some other 
methods, namely shift-add through all segments of data (see, Barret et 
al., 2005), provide more efficient analysis of the frequency correlation 
but do not keep the information about the probability of the detection 
(i.e., the frequency distribution). 



of a test particle is not closed and the frequencies of the az- 
imuthal, radial and vertical "quasielliptic" orbital motion differ. 
For a given axially symetric spacetime with the line element 

ds 2 = g„ dt 2 + g rr dr 2 + g m dQ 2 + dcf> 2 + g M d<f> dt + 



(1) 



the relevant angular velocities read (e.g., Aliev & Galtsov, 1981; 
Okazaki et al., 1987; Stella & Vietri, 1998; Abramowicz et al., 
2003a) 3 



Ok 



2 _ igu + Q ± g«pf ( d 2 U 



where i e (r, 8) and U is an effective potential 
U(r, 0,€):= g" - 21 g"* + t 2 g'f'f, 



(2) 
(3) 

(4) 



with { denoting the specific angular momentum of the orbiting 
test particle 



£ = -ua,Iu, . 



(5) 



In the following we consider Keplerian motion and I = 
lic(r, 0). The motion is then described by the Keplerian frequency 
v K and radial and vertical epicyclic frequencies v r , v e . 

Due to the inequality between the azimuthal and radial fre- 
quency, the eccentric orbits waltz at the periastron precession 
frequency (e.g., Misner et al., 1973) 

v P = v K - v T , (6) 

and in addition the orbits tilted relative to the equatorial plane 
of the spinning central mass wobble at the nodal (often called 
Lense-Thirring) precession frequency (e.g., Misner et al., 1973) 



ve- 



il) 



The periastron precession frequency (6) corresponds to the 
period in which the test particle "quasielliptic" trajectory peri- 
astron oscillates and the nodal (Lense-Thirring) precession fre- 
quency (7) corresponds to the period in which the declination 



3 In the text above we mark the frequencies expressed in standard (SI) 
physical units by asterisk. Hencefort we already use the geometrical 



units (c = 1 = G; M = GM'/c 1 



t = cf). 
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of the quasiellipse plane oscillates. Both the declination of the 
quasiellipse plane and position of the periastron then reach the 
initial state simultaneously in the period characterized by the fre- 
quency 



Vg 



(8) 



Therefore, for the purposes of our paper, we call this frequency 
total precession frequency. 4 

2.1. The Hartle-Thorne metric 

In this paper we describe the external neutron star spacetime 
using the Hartle-Thorne metric (Hartle & Thorne, 1968) which 
represents the exact solution of vacuum Einstein field equations 
for the exterior of rigidly and relatively slowly rotating, station- 
ary and axially symmetric body. The metric is given with ac- 
curacy up to the second order terms in the body's dimensionless 
angular momentum j = J/M 2 , and first order in its dimensionless 
quadrupole moment q = -Q/M 3 . The explicit form of formulae 
(2) and (3) derived by Abramowicz et al. (2003a), which we use 
in a slightly modified form, is given in Appendix A. 

3. Testing the multiresonant hypothesis 

3.1. Frequency identification 

Usually the n : m orbital resonant models considering a non- 
linear resonance between Keplerian and/or epicyclic frequencies 
(see, e.g. Abramowicz et al., 2004) identify the resonant eigen- 
frequencies vf , as 



v£ = v r (r n:m ), v° = v v (r n:m ), v v e [v g , v K ] 



(9) 



where n, m are small natural numbers and r n:m is the radius fixed 
by the condition 



Vy(r„:m) _ «_ 
Vr(r„ : „,) m 



(10) 



In the case of a considerably weak forced or parametric non- 
linear resonance (Landau & Lifshitz, 1976), the upper and lower 
observed QPO frequencies v L and v v are related to the resonant 
eigenfrequencies either directly 



. o . o 

Vr = V, , V„ = V„, 



or as their linear combinations 

v L = av Q L , v„=/3v", 



(11) 



(12) 



where a and ft are small integral numbers. This property was 
utilized to estimate the spin of microquasars displaying con- 
stant twin peak QPOs from resonant models (Abramowicz & 
Kluzniak, 2001; Torok et al., 2005). 

In general case of a system in a non-linear resonance, the ob- 
served frequencies differ from resonance eigenfrequencies by a 
frequency corrections proportional to the square of small dimen- 
sionless amplitudes (Landau & Lifshitz, 1976). It was shown 
(Abramowicz et al., 2003b; Rebusco, 2004) that a resonance 
characterized by one pair of eigenfrequencies may reproduce the 
whole range of frequencies observed in a neutron star source. 
Later Abramowicz et al. (2005) considered the idea of one 



4 Note that in general the repeating of initial periastron position and 
orbit declination do not guarantee repeating of the initial test particle 
position at the same time. 



eigenfrequency pair (so called resonant point in the frequency- 
frequency plane) common for a set of neutron star sources. They 
shown that for weakly coupled non-linear resonance the upper 
and lower frequency observed in a source should be linearly 
correlated. They also found that coefficients of linear fits well 
approximating individual sources are antic orrelated which was 
in a good accord to the theory they presented and justified the 
hypothesis of one eigenfrequency-pair. On the other hand this 
approach, incorporating certain difficulties (e.g., the extremely 
large extension of the observed frequency range), is not proved 
yet, and some observational facts like the multipeaked ratio dis- 
tribution suggest that more then one resonant points may be re- 
sponsible for the almost linear observed frequency correlation. 

In next we focus on the hypothesis of more resonant points 
corresponding to different instances of one resonance and sup- 
pose that the observed frequencies are close to the resonance 
eigenfrequencies, i.e. that the observed frequency correlation 
follows the generic relation between resonant eigenfrequencies, 

v L ~ v°, v v ~ v°. (13) 

For the whole applicable range of the internal angular mo- 
mentum j of the Hartle-Thorne spacetimes we checked that the 
ratio between the Keplerian (or vertical epicyclic) frequency and 
radial epicyclic frequency monotonically increases with decreas- 
ing radius r whereas the Keplerian (vertical epicyclic) frequency 
increases (see also Torok & Stuchlfk, 2005). 

In other words, for the models (9) considering resonance 
between Keplerian (vertical epicyclic) frequency and radial 
epicyclic frequency satisfying relation (13), the ratio of observed 
frequencies should increase with increasing QPO frequency, but 
that is opposite to what is observed. 

However, the relations (10-12) are not the only possible in 
the framework of resonance models. Bursa (2005) discussed so 
called vertical precession resonance introduced in order to match 
the spin estimated from fits of the X-ray spectral continua for 
the microquasar GRO Jl 655-40. The resonance should occur be- 
tween the vertical epicyclic frequency and the periastron preces- 
sion frequency fulfilling the relation 

v°(r) = v P (r) = v K (r) - v r (r), v°(r) = v e (r), (14) 

for a particular choice of the resonant radius r defined by the 
condition v u = 3/2vi. 

As noticed in Torok et al. (2007), for the Schwarzschild 
spacetime the relations (14) coincide with those following from 
the relativistic precession model: 

v° L (r) = v P (r) = v K (r) - v r (r), v°(r) = v K (r). (15) 

Opposite to the relations (9) the two relationships (14,15) as well 
as the other two relationships 



v» = v fl (r)-v r (r), 



y» = v fl (r), 



v?(r) = v T (r) = vg(r) - v r , v„°(r) = v K (r) 



(16) 
(17) 



imply the increase of for increasing v'\ 

In the following subsection we fit the QPO frequencies ob- 
served in 4U 1636-53 by the four different frequency relation- 
ship (14-17), testing the hypothesis that an appropriate reso- 
nance may be responsible for all the observed datapoints. 

3.2. Individual fits 

The fits are realized in two steps. In order to obtain a rough 
scan we calculated frequency relations (14-17) in the Hartle- 
Thorne metric for the range of the mass M e 1-3Mq, the inter- 
nal angular momentum j e 0-0.5 and a physically meaningful 
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Fig. 2. a) The inverse quality measure^ 2 of the datapoints fits by relation (15) -"Stella" as a function of the neutron star mass M and the angular 
momentum j. The displayed individual values of x 2 correspond to the quadrupole momentum q giving the best fit for the fixed pairs (M, j). 
b) The same for the relation (14)-"Bursa". c) The same but for the relation (17)- total precession. For all the three relations, the^ 2 minima form 
"valleys" with rather sharp roots, d) Profiles of the individual valleys. The cross-hair marks denote the value of the global minima. For relations 
"Stella" and "Bursa" the value of x 2 local minimum at the bottom of the valley decreases with the increasing angular momentum j. On the other 
hand for the total precession relation it increases with the increasing angular momentum j. 



quadrupole momentum q with a step equivalent to the thousand 
points in all three quantities, i.e., four 3-dimensional maps each 
having 10 points. Then, for each pair (M,j), we keep the value 
of the quadrupole momentum q which gives the lowest x 2 with 
respect to the observed datapoints. For the Schwarzschild space- 
time (q — j = 0), when relations (14-17) merge, the best fit is 
reached for the mass M = 1 .11 Mq, with a^ 2 = 400 ~ 20 d.o.f. 
Because for rotating configurations the relation (16) gives rather 
less interesting results, we do not consider this relation in the 
following. 

Figures 2 a, b, and c show the lowest^ 2 associated with rela- 
tions (15) -"Stella", (14) -"Bursa", and (17)- total precession, 
for each pair (M, j). Having a rough clue given by Figures 2a,b, 
and c we searched for local x 1 minima using the Marquardt- 
Levenberg non-linear least squares method (Marquardt, 1963). 
Particular minima we have found are denoted in Figure 2d. 

The detailed analysis shown that both the relations (14) and 
(15) match the observational data most likely for relatively high 
angular momentum close to j ~ 0.5 and the central mass M ~ 2.4- 
2.8Mq ("Stella"), and M-2.4M© ("Bursa") respectively, reach- 
ing the value of^ 2 ~ 15 d.o.f., whereas the quadrupole momen- 
tum is rather close to the Kerr value q — a 1 — j 1 (q~0.23, 0.25). 

Relation (17, total precession) gives the best fit with the re- 
markable x 1 ~ 3 d.o.f. for angular momentum j ~ 0.05 and cen- 



tral M ~ 1 .76 Mq, again with the quadrupole momentum close 
to the Kerr value (^~0.029) 5 . For this relationship the quality of 
the fit is then very close to x 2 ~ 3 d.o.f. in rather large range of 
the central mass M~ 1.76-1.84Mq and the angular momentum 
j~ 0.05-0.1. 

4. Discussion and conclusions 

Figures 3a,b, and c show the best fits reached by relations (14- 
17). The relevant best fits properties are summarized in the Table 
1. Apparently, the identification of the lower kHz QPOs with 
total precession frequency and the upper kHz QPOs with the 
Keplerian frequency 6 

(r) - v T (r) = v e (r) - v r (r), v°(r) = v K (r) 

reaches a substantially better quality of the fit (x 2 = 3 d.o.f.) 
then the other three possibilities. The corresponding value of 



5 Our results are thus compatible with expectation that gravitational 
field of settled down neutron stars can be well described by quasi-Kerr 
spacetime. 

6 As previously stressed, our paper focuses on the hypothesis of one 
resonance which occurs at different resonant points. The more general 
multiresonant idea, including the possibility of several resonances shar- 
ing the common radius, relevant for compact objects with special values 
of the angular momentum, is discussed in Stuchlfk et al. (2007). 
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Fig. 3. a) The best fit by relation (15) -"Stella", b) The best fit by relation (14)-"Bursa". c) The best fit by relation (17)-total precession, d) 
Adopted from Torok & Barret (2007). The difference of the lower and upper QPOs rms amplitudes plotted as a function of the lower QPOs 
frequency (black points). Assuming zero amplitude of the missing QPOs mode, the quantity Ar is for single peaks (grey points) given by an 
absolute value of their rms amplitude. To plot Ar corresponding to the upper single QPO we used the linear frequency-frequency fit. Dotted 
vertical lines denote position of the qualitative change in the rms amplitude interrelationship. 



the angular momentum j ~ 0.05-0.1 and mass M ~ 1.8 are 
not excluded by any study published so far (Strohmayer & 
Markwardt, 2001; Giles et al., 2002) and this frequency identi- 
fication appears more likely than the relationship corresponding 
to the relativistic precession model of Stella & Vietri (M ~ 2.4- 
2.8M©, 7-0.4-0.5). 

The total precession frequency v T given by (8) corresponds 
to the period in which the declination of the free test particle 
quasiellipse plane and the periastron reach simultaneously the 
initial state. Consideration of a hot spot with characteristic fre- 
quencies (17) may represent a kinematic QPOs model very close 
(but obviously not identical) to those of Stella & Vietri. 7 

However, the observed ratio clustering and rms amplitude 
difference behaviour suggest the existence of a resonance be- 
tween the lower and upper QPOs modes. Notice also that the 
position of the 3/2 and 5/4 resonant points implied by the total 
precession relationship (i.e., intersections of the best datapoints 
fit with reference 3/2 and 5/4 lines) well coincide with frequen- 
cies given by the change of the rms-amplitude difference sign 



(see Figure 3 and Table 1). Therefore, our results indicate that 
the resonance may occur between the Keplerian frequency of 
the trajectory and the total precession frequency corresponding 
to the periodicity of the trajectory shape. 8 The concrete physical 
mechanism of such a resonance remains to be the subject of a 
future and rather larger research since, beyond the hotspot inter- 
pretation, the Keplerian and total precession frequency may also 
correspond to some disk oscillation modes. 

The results we have obtained so far indicate that Keplerian 
and total precession frequency (17) match the kHz QPOs fre- 
quencies very well at least in several atoll sources (Bakala et al., 
2007, in preparation). 

On the other hand, the applicability of the relation (17) to 
the Galactic microquasar sources poses an open question be- 
cause it implies the central black hole spin similar to those of 
Stella & Vietri (a ~ 0.3-0.5), which in the case of microquasar 
GRS 1915+105 contradicts the recent results of a fitting the 
spectral continua (a ~ 0.7-1, McClintock et al., 2006; Middleton 
et al., 2006). 



7 One may also argue that if the light curve is somehow simultane- 
ously modulated by the periastron and Lense-Thirring precession, the 
beat frequency should appear in its power spectra. 



8 For the perfect free particle motion, if the Keplerian and total pre- 
cession frequency form rational fractions, the trajectory is selfrepeating 
(i.e., closed). 
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Table 1. The properties of best fits depicted in Figure 3a,b, and c (in the case of the totall precession we display also properties of the fit for the 
angular momentum j = 0.1). The uncertainities in the fit parameters are the standard^ 2 + 1 errors. Quantity Aq characterizes the deviation of the 
quadrupole momentum from the value corresponding to the Kerr spacetime. Quantities Av L characterize the deviation of best fit intersections with 
reference 3/2 and 5/4 lines from the value given by the qualitative change in the rms amplitudes behaviour. 
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Appendix A: Formulae for orbital geodesic frequencies in the Hartle-Thorne metric 

After Abramowicz et al. (2003a), the Keplerian orbital angular velocity and the radial and vertical epicyclic angular velocities can 
be expressed in terms of the Hartle-Thorne metric parameters M, j, , q in the following form. 

The angular velocity for corotatingcircular particle orbits reads 



where 



Mi/2 

-3/2 



M 3/2 

1- j—+f + qF 2 n (r) 



F^(r) = [48 M 1 - 80 M 6 r + 4 M 5 r - 18 M 4 r 3 + 40 M 3 r 4 + 10 M 2 r 5 
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The epicyclic frequencies of circular geodesic motion are given by formulae 
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For completness, the relation determining the marginally stable circular geodesic reads 
r ms = 6 M 
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